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Two models for energy and charge transport and storage in biomolecules are considered. A 
model based on the discrete nonlinear Schrodinger equation with long-range dispersive interactions 
(LRI's) between base pairs of DNA is offered for the description of nonlinear dynamics of the DNA 
molecule. We show that LRI's are responsible for the existence of an interval of bistability where 
two stable stationary states, a narrow, pinned state and a broad, mobile state, coexist at each value 
of the total energy. The possibility of controlled switching between pinned and mobile states is 
demonstrated. The mechanism could be important for controlling energy storage and transport in 
DNA molecules. Another model is offered for the description of nonlinear excitations in proteins 
and other anharmonic biomolecules. We show that in the highly anharmonic systems a bound state 
of Davydov and Boussinesq solitons can exist. 
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I. INTRODUCTION 

Understanding how biological macromolecules (pro- 
teins, DNA, RNA, etc.) function in the living cells re- 
mains the major challenge in molecular biology. One 
of the most important questions is the mechanism of 
gene expression. The expression of a given gene involves 
two steps: transcription and translation. The transcrip- 
tion includes copying the linear genetic information into 
the messenger ribonucleic acid (mRNA). The informa- 
tion stored in mRNA is transferred into a sequence of 
aminoacids using the genetic code. mRNA is produced 
by the enzyme RNA-polymerase (RNAP) which binds 
to the promoter segment of DNA. As a result of the in- 
teraction between RNAP and promoter of DNA the so- 
called "bubble" (i.e. a state in which 10-20 base pairs 
are disrupted) is formed. The disruption of 20 base pairs 
corresponds to investing some 100 kcal/mole (0.43 eV) 

0. 

In the framework of a linear model the large-amplitude 
motion of the bases was supposed to occur due to an in- 
terference mechanism [||. According to this model ener- 
getic solvent molecules kick DNA and create elastic waves 
therein. As a result of the interference of two counter 



propagating elastic waves, the base displacements may 
exceed the elasticity threshold such that DNA undergoes 
a transition to a kink form which is more flexible. A sim- 
ilar approach was also proposed in Refs. j|,§]|- The linear 
elastic waves in DNA are assumed to be strong enough to 
break a hydrogen bond and thereby facilitate the disrup- 
tion of base pairs. In spite of the attractiveness of this 
theory, which gives at least a qualitative interpretation 
of the experimental data |J , there are the following fun- 
damental difficulties which to our opinion are inherent 
in the linear model of the DNA dynamics: (i) The dis- 
persive properties (the dependence of the group velocity 
on the wave length) of the vibrational degrees of free- 
dom in DNA will cause spreading of the wave packets 
and therefore smear the interference pattern. Further- 
more, it has been shown [0 that the amplitudes of the 
sugar and the base vibrations are rather large even in a 
crystalline phase of DNA. Since the large-amplitude vi- 
brations in the molecules and the molecular complexes 
are usually highly anharmonic their nonlinear properties 
can not be ignored; (ii) Molecules and ions which exist in 
the solution permanently interact with DNA. These in- 
teractions are usually considered as white noise and their 
influence is modelled by introducing Langevin stochas- 



tic forces into the equations describing the intramolec- 
ular motion. It is well known g] that stochastic forces 
provide relaxation of linear excitations and destroy their 
coherent properties. Equivalently the coherence length 
(the length of the concerted motions) rapidly decreases 
with increasing temperature; (iii) DNA is a complex sys- 
tem which has many nearly isoenergetic ground states 
and may therefore be considered as a fluctuating aperi- 
odic system. DNA may have physical characteristics in 
common with quasi-one-dimensional disordered crystals 
or glasses. However, it is known [§) that the transmission 
coefficient for a linear wave propagating in a disordered 
chain decreases exponentially with the growth of the dis- 
tance (Anderson localization). In this way it is difficult 
to explain in the framework of linear theory such a phe- 
nomenon as an action at distance, where concerted mo- 
tion initiated at one end of a biological molecule can be 
transmitted to its other end. 

The above mentioned fundamental problems can be 
overcome in the framework of nonlinear models of DNA. 
Nonlinear interactions can give rise to very stable excita- 
tions, called solitons, which can travel without changing 
their shape. These excitations are very robust and impor- 
tant in the coherent transfer of energy M . For realistic 
interatomic potentials the solitary waves are compressive 
and supersonic. They propagate without energy loss, and 
their collisions are almost elastic. 

Nonlinear interactions between atoms in DNA can 
give rise to intrinsically localized breather-like vibra- 
tion modes |10| , pT[ . Such localized modes, being large- 
amplitude vibrations of a few (2 or 3) particles, can facil- 
itate the disruption of base pairs and in this way initiate 
conformational transitions in DNA. These modes can oc- 
cur as a result of modulational instability of continuum- 
like nonlinear modes p2| , which is created by energy ex- 
change mechanisms between the nonlinear excitations. 
The latter favors the growth of the large excitations [[13| . 

Nonlinear solitary excitations can maintain their over- 
all shape on long time scales even in the presence of the 
thermal fluctuations. Their robust character under the 
influence of white noise was demonstrated |14| and a sim- 
plified model of double-stranded DNA was proposed and 
explored. Quite recently the stability of highly localized, 
breather-like, excitations in discrete nonlinear lattices 
under the influence of thermal fluctuations was investi- 
gated pi . It was shown that the lifetime of a breather 
increases with increasing nonlinearity, and in this way 
these intrinsically localized modes may provide an exci- 
tation energy storage even at room temperatures where 
the environment is strongly fluctuating. 

Several theoretical models have been proposed in the 
study of the nonlinear dynamics and statistical mechan- 
ics of DNA (see the very comprehensive review JT(|). A 
particularly fruitful model was proposed by Peyrard and 
Bishop [n7[ and Techera, Daemen and Prohofsky ffq] . In 
the framework of this model the DNA molecule is con- 
sidered to consist of two chains that are transversely cou- 
pled. Each chain models one of the two polynucleotide 



strands of the DNA molecule. A base is considered 
to be a rigid body connected with its opposite partner 
through the hydrogen-bond potential V{u n ), where u n is 
the stretching of the bond connecting the bases, n, and 
n = 0,±I,±2, .. is labelling the base-pairs. The stretch- 
ing of the n'th base-pair is coupled with the stretch- 
ing of the m'th base-pair through a dispersive potential 
U(u n , u m ). The process of DNA denaturation was stud- 
ied J17| , [19[ under the assumption that the coupling be- 
tween neighboring base-pairs is harmonic U(u n ,u n +i). 
An entropy-driven denaturation was investigated |2C| ] 
taking into account a nonlinear potential U(u n , w n +i) be- 
tween neighboring base-pairs. The Morse potential was 
chosen |I^,|I^,^] as the on-site potential V(u n ), but also 
the longitudinal wave propagation and the denaturation 
of DNA has been investigated [H using the Lennard- 
Jones potential to describe the hydrogen bonds. 

In the main part of the previous studies the disper- 
sive interaction U was assumed to be short-ranged and 
a nearest-neighbor approximation was used. It is worth 
noticing, however, that one of two hydrogen bonds which 
is responsible for the interbase coupling, namely the hy- 
drogen bond in the N — H...0 group, is characterized by a 
finite dipole moment. Therefore, a stretching of the base- 
pair will cause a change of the dipole moment, so that the 
excitation transfer in the molecule will be due to transi- 
tion dipolc-dipole interaction with a 1/r 3 dependence on 
the distance, r. It is also well known that nucleotides in 
DNA are connected by hydrogen-bonded water filaments 
f2lj , f22[ . In this case an effective long-range excitation 
transfer may occur due to the nucleotide-water coupling. 

In the last few years the importance of the effect of 
long-range interactions (LRI) on the properties of non- 
linear excitations was demonstrated in several different 
areas of physics [23]-|2^|. Quite recently p9[ , we proposed 
a new nonlocal discrete NLS model with a power depen- 
dence on the distance of the matrix element of the disper- 
sive interaction. It was found that there is an interval of 
bistability in the NLS models with a long-range disper- 
sive interaction. One of these states is a continuum-like 
soliton and the other is an intrinsically localized state. 

Another model, which we consider in this paper, de- 
scribes an excitation energy (or excess electron) trans- 
port in a highly anharmonic molecular chain. This line 
of investigation is rather new, since until recently most 
attention has been concentrated on strongly elastic sys- 
tems, such as the excitons in the a-helix proteins. The 
subsonic Davydov solitons can exist in this case, and the 
harmonic approximation is enough for a consideration 
of the lattice deformation caused by the exciton Q.|fJ . 
But if the lattice is not supposed to be rigid or if the mo- 
tion of the excitation at supersonic velocities is studied, 
the deformation cannot be considered small and anhar- 
monic terms must be taken into account. For example, 
the anharmonic interactions play an important role in 
the dynamics of DNA. It is possible that they might 
concentrate vibrational energy in DNA into soliton-likc 
objects MM. 



The anharmonic Davydov model of the molecular 
chain has already been much studied during the last fif- 
teen years |3(J,|33| |39|. It is well known that there exist 
solitons of Davydov type and of Boussinesq type in the 
system. Perez and Theodorakopoulos E3] investigated 
numerically an interaction between these solitons in the 
a-helix proteins, and concluded that it inhibits the sta- 
bility of the Davydov vibrational soliton at room temper- 
ature. But the effective anharmonicity parameter in that 
system is rather small. The recent investigations on the 
subject pq-p9] demonstrated that in highly anharmonic 
chains (e.g. in the problem of excess electron transfer in 
a-helix proteins), the interaction between Davydov and 
Boussinesq solitons changes its sign, and the solitons can 
form a bound state. 

The outline of the paper is the following. In Sec. II 
we investigate the effects of long-range interactions on 
the nonlinear dynamics of the two-strand model of DNA. 
First, we reduce the equations of motion of the model to 
the nonlocal discrete nonlinear Schrodinger (NLS) equa- 
tion. Then, we present the analytical theory and the re- 
sults of numerical simulations of the discrete NLS model 
with a long-range dispersive interaction. We obtain sta- 
tionary localized solutions, demonstrate the existence of 
the bistability phenomenon, and discuss their stability. 
The investigation of the switching between bistable states 
completes the section. We show that a controlled switch- 
ing between narrow, pinned states and broad, mobile 
states with only small radiative losses is possible when 
the stationary states possess an internal breathing mode. 
In Sec. Ill we investigate the anharmonic Davydov model 
of the molecular chain and show that it reduces to the 
Henon - Heiles system Ell], which is completely inte- 
grable at some fixed value of anharmonicity. There are 
three types of stationary solitons in this case: Boussi- 
nesq, Davydov, and a new type, namely a supersonic 
two-bell shaped soliton. We calculate the energies of the 
solitons. Then, we develop a variational approach for the 
non-integrable case. We show that in the case of highly 
anharmonic chains, the two-bell shaped soliton is an os- 
cillating bound state of Davydov and Boussinesq solitons. 
This bound state is caused by the excitation (or electron) 
tunnelling in the effective two-well potential which is cre- 
ated by the exciton (electron) - phonon interaction and 
anharmonic terms in the lattice potential. On the con- 
trary, at small anharmonicity Davydov and Boussinesq 
solitons repel each other, and the bound state cannot ex- 
ist. Finally, we test the dynamical stability of the bound 
state at various values of anharmonicity numerically. 



II. EFFECTS OF LONG-RANGE DISPERSION: 
BISTABILITY OF LOCALIZED EXCITATIONS 



A. System and equations of motion 

In this section we study the two-strand model of DNA 
which is described by the Lagrangian 



L = T-U B p- U lr 



where 



T = 



M 



^r 



du r 



(1) 



(2) 



2 ^ V dt 

n 

is the kinetic energy and M is the mass of the base-pair, 

u BP =j2 v M ( 3 ) 

n 

is the potential energy which describes an intrabase-pair 
interaction, and 

1 



U L r = 



sr sr 



4 Z^ Z^ 



Jn,m\U"m ^n) 



(4) 



is the long-range dispersive interbase-pair interaction of 
the stretchings. In Eqs. (H)-(Q) n and m are base-pair 
indices. The value u n = of the base-pair stretching 
corresponds to the minimum of the intrabase-pair poten- 
tial V(u n ). We investigate the model with the following 
power dependence on the distance (assuming that the 
lattice constant equals unity) of the matrix element of 
the base elastic coupling 

Jn,m = -, jT , (5) 

1 7i — m\ s 

where the constant J characterizes the strength of 
the coupling, and s is a parameter being introduced 
to cover different physical situations including the 
nearest- neighbor approximation (s = oo), quadrupole- 
quadrupole (s = 5) and dipole-dipole (s = 3) interac- 
tions. In the case of the DNA molecule the long-range in- 
teraction is evolved from the existence of charged groups 
in the DNA molecule. To take into account a possibil- 
ity of screening of the interaction or an indirect coupling 
between base-pairs (e.g. via water filaments), we shall 
consider also the case when the matrix element of the 
base elastic coupling has the Kac-Bakcr [42(43] form 



Jn.m — " G 



-/3|n— m| 



where (3 is the inverse radius of the interaction. 
Assuming that 



d 2 V(u n ) 



dul 



» 



d' J V(u n ) 



«„=o 



dui 



(6) 



(') 



u„=0 



for j = 3, 4, ..., that is the anharmonicity of the intrabase- 
pair potential is rather small, we shall use a rotating-wave 
approximation 

-iflt 



il> n e~ 



+ c.c. 



(8) 



M„=0 



is the frequency of the 



where il = ^/V"{u n )/M 

harmonic oscillations and tp n (t) is the complex amplitude 



which is supposed to vary slowly with time. Inserting Eq. 
(|J) into Eqs. (§])-([|) and averaging with respect to the 
fast oscillations of the frequency fi, we conclude that the 
effective Lagrangian of the system can be represented in 
the form 



^=lJ2(^n-^n)-H 



(9) 



where the dot denotes the differentiation with respect to 
the rescaled time r = t/(2MCl). Here 

H=U BP +U LR (10) 

is the effective Hamiltonian of the system, where 



n,m(n^rn) 

is the effective dispersive energy and 

n r 27r/n 



(ii) 



•*-?&£ 



V(ip n e- int + c.c)dt 

Mn 2 \ip n \ 2 ) (12) 



is the effective intrabase-pair potential. Usually either a 
Morse potential [|l7|,[l9| or a Lennard- Jones potential |1J] 
is used to model the hydrogen bonds. With these poten- 
tials however it is very complicated to obtain any analyt- 
ical results. Therefore, to gain insight into the problem 
we shall use a simplified nonlinear potential in the form 



Ubp — — 



(a + 1)^ 



IV'nl 



(13) 



where the degree of nonlinearity a is a parameter which 
we include to have the possibility to tune the nonlinearity 
as well. 

From the Hamiltonian (fill) we obtain the equation of 
f¥H — 

motion iip n — - — for the wave function tynij) m f ne 



form 



dip* n 



#n+ Yl J n,m(^m-^n) + \^n\ 2a ^n^0- (14) 



i(m^r. 



The Hamiltonian TC and the number of excitations 



^=Ew 



(15) 



are conserved quantities. 

Thus we have reduced the equations of motion of the 
initial model of DNA dynamics to the nonlocal discrete 
nonlinear Schrodinger (NLS) equation studied quite re- 
cently in detail paH. In what follows we shall review 
the obtained results in application to the DNA molecule. 



Looking for the stationary solutions of Eq. ( |l4| ) of the 
form 



ipn = <t>n exp(iAr) 



(16) 



with a real shape function <p n and frequency A, we obtain 
the equation for <p n in the form 



A0„ = 22 Jn,m(4>m 



t>n) + C 



2(7 + 1) 



(17) 



with J n , m described by Eqs. (|) or (§). Thus Eq. @ is 
the Euler-Lagrange equation for the problem of extrem- 
izing TL under the constraint ./V = constant. 



B. Bistability of stationary states 

To develop a variational approach to the problem we 
use an ansatz for a localized state in the form 



(f> n = v7Vtanhaexp(— a\n\) , 



(18) 



where a is a trial parameter. This ansatz is chosen to 

satisfy automatically the normalization condition (|15| ) , so 

that the problem of extrcmizing TL under the constraint 

N = constant is reduced to the problem of solving the 

• dH n 
equation — — = 0. 
da 
Inserting the trial function (|18j) into the Hamiltonian 

given by Eqs. (|10|), (|ll|), ( |l^ ) and (||) and evaluating the 

discrete sums which enter in these equations (see p9[ | for 

details), we get the dispersive part of the Hamiltonian 



U LR = 2NJ{((s) - tanh(a)F(e~", s - 1) 
-F(c~ a ,s)} 

and the intrabase-pair potential 

N a+1 



(19) 



Ubp — 



u, 



(7+1 

with f a = tanh ff+1 (a) coth[(cT + l)a] 



(20) 



where 

oo 



and F(z,s) = J2( zn / nS ) ( 21 ) 



n=l 



are Riemann's zeta function and Jonqiere's function, re- 
spectively. 

According to the variational principle we should satisfy 

the condition — — = which yields 



dot 

N a = 2(cr + 1) J (tanh(a)F(e- Q , s - 2) 

+ tanh 2 (a)F(e- Q ,s-l)) 



(22) 



dU 

da 



As a direct consequence of Eq. (|l7J), the frequency A can 
be expressed as 



A = --{U LR + 2U B p) 



(23) 



with Ulr and Ubp being defined by Eqs. ( Jiq| ) and (pO|) . 

The results of the above developed variational ap- 
proach are plotted in Figure |l|, which shows the energy 
of the localized excitations as a function of N and s for 
the particular case a = 1. 




FIG. 1. The dependence of the energy H of stationary soli- 
ton solutions versus the excitation number N and the disper- 
sion parameter s, obtained from the variational approach for 

C7=l. 



One can see that there is a critical value s cr = 2.72 
of the dispersion parameter s which separates two qual- 
itatively different regions. The region s > s cr is char- 
acterized by the usual monotonic behavior of the soliton 
energy vs number of excitations N . Thus the main fea- 
tures of all discrete NLS models with dispersive inter- 
action J ntm decreasing faster than \n — 77i| _Scr coincide 
qualitatively with the features obtained in the nearest- 
neighbor approximation where only one stationary state 
exists for any N. But for 2 < s < s cr the bifurcation 
"swallow tail" takes place introducing a multistability 
into the system. One can see that in this case there is 
an energy interval where for each value of Ti there exist 
three stationary states with different excitation numbers. 
The direct numerical solution of Eq. (|17| ) validates this 
conclusion and gives the precise critical value s cr = 3.03. 
It is noteworthy that similar results are obtained for the 
Kac-Baker dispersive interaction (13). In this case the 
multistability takes place for ft < 1.70. 



5 
4 

N 3 
2 
1 




\ ' ■ ' 


- \ 




<>V* ,.i— ^**^"^ 


I / jf 


1/ 


1 1 1 



0.0 



0.5 



1.0 
A 



1.5 



2.0 



FIG. 2. Number of excitations, TV, versus frequency, A, nu- 
merically from Eq. ( |l7| ) for a = 1 and s = oo (full), 3 (dotted), 
2.3 (dashed), 2 (short- long-dashed), and 1.9 (long-dashed). 
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FIG. 3. Shapes of the three stationary states for s = 2.3 
and N = 3.2. The stable: A = 0.11 and A = 0.98 (full). The 
unstable: A = 0.43 (dashed). 

Figure shows the dependence N(A) obtained numer- 
ically for a = 1. This dependence demonstrates the same 
transition from the monotonic behavior for s > s cr to the 
non-monotonic (TV-type) behavior for s < s cr and assists 
to determine the stability of the stationary states. Ac- 
cording to the theorem which was recently proven ||45| , 
the necessary and sufficient linear stability criterion for 
the stationary states is 



dN d v-^ 9 



(24) 



Therefore we can conclude that only two stationary 
states are stable. The third state, with intermediate 
value of A and maximal value of energy, is unstable. Thus 
it is reasonably to talk about bistability of the system. 

The observed bistability is very similar to the one re- 
cently observed p5 46|, where the nearest-neighbor case 
with an arbitrary degree of nonlinearity a was studied. 



The bistability appears in this case for a above a certain 
critical value. 

Figure g shows that the shapes of these solutions dif- 
fer significantly. The low-frequency states are wide and 
continuum-like, while the high-frequency solutions repre- 
sent intrinsically localized states with a width of a few 
lattice spacings. 

Now we turn to discuss stationary states of the discrete 
NLS model given by Eq. ( |17|) with arbitrary degree of 
nonlinearity. The main properties of the system remain 
unchanged, but the critical value s cr of the dispersion pa- 
rameter is now a function of a. The results of analytical 
consideration confirmed by simulation show that s cr in- 
creases with increasing a. In particular, for a > 1.4 (the 
value at which the discrete symmetrical ground state can 
be unstable in the nearest-neighbor approximation |45[ ) 
the bistability in the nonlinear energy spectrum occurs 
even for s < 6. 



C. Switching between bistable states 

Having established the existence of bistable station- 
ary states in the nonlocal discrete NLS system, a natural 
question that arises concerns the role of these states in 
the full dynamics of the model. In particular, it is of in- 
terest to investigate the possibility of switching between 
the stable states under the influence of external pertur- 
bations, and what type of perturbations that could be 
used to control the switching. Switching of this type is 
important in the description of nonlinear transport and 
storage of energy in biomolecules like the DNA, since a 
mobile continuum-like excitation can provide action at 
distance, while the switching to a discrete, pinned state 
can facilitate the structural changes of the DNA |||l6) . 

An illustration of how the presence of an internal 
breathing mode can affect the dynamics of a slightly per- 
turbed stable stationary state is given in Figures [| and ||. 
To excite the breathing mode we apply a spatially sym- 
metrical, localized perturbation, which we choose to con- 
serve the number of excitations in order not to change the 
effective nonlinearity of the system. The simplest choice, 
which we have used in the simulations shown here, is to 
kick the central site no of the system at t = by adding 
a parametric force term of the form 0S n ^ no 5(T)ip n (T) to 
the left-hand-side of Eq. ( |l4| ) . All details and a physical 
motivation of the appearance of such kind of parametric 
kick are described in Ref. p4[ . As can be easily shown, 
this perturbation affects only the site Uq at r = 0, and re- 
sults in a 'twist' of the stationary state at this state with 
an angle 8, i.e. tp no (0) = </>„„ e l6 '. The immediate conse- 
quence of this kick is, as can been deduced from the form 

of Eq. (|1J), that — (|^ no | 2 ) will be positive (negative) 

when 8 > ( 8 < 0). Thus, to obtain switching from 
the continuum-like state to the discrete state we choose 
8 > 0, while we choose 8 < when investigating switch- 
ing in the opposite direction. We find that in a large part 



of the multistability regime there is a well-defined thresh- 
old value 8th, such that when the initial phase torsion is 
smaller than 8 t h, periodic, slowly decaying 'breather' os- 
cillations around the initial state will occur, while for 
strong enough kicks (phase torsions larger than 8th) the 
state switches into the other stable stationary state. 
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FIG. 4. Switching from continuum-like to discrete state 
for the NLS model with the Kac-Baker long-range interac- 
tions (H) where (3=1. The initial state cf> n has the fre- 
quency A ~ 0.310 and N = 3.6. We plot a time evolution 
of | - i/ , n ('r)| 2 when a phase torsion is applied to the center 
site with 9 = 0.261 (lower curve) and 9 = 0.262 (upper 
curve), respectively; inset shows time evolution of IV'nMI 2 
for 9 = 0.262. 
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FIG. 5. Switching from discrete to continuum-like state 
for the NLS model with the Kac-Baker long-range interac- 
tions (H) where /3 = 1. The initial state <j> n has the fre- 
quency A ~ 1.423 and A" = 3.6. We plot a time evolution 



Of \i>n 



when a phase torsion is applied to the center 



site with 9 — —0.490 (upper curve) and 9 = —0.491 (lower 
curve), respectively; inset shows time evolution of |V'»i(' r )| 2 f° r 
9 = -0.491. 

It is worth remarking that the particular choice of per- 
turbation is not important for the qualitative features 
of the switching, as long as there is a substantial over- 
lap between the perturbation and the internal breathing 
mode. We also believe that the mechanism for switching 
described here can be applied for any multistable system 
where the instability is connected with a breathing mode. 



For example, we observed [J47| a similar switching behav- 
ior in the nearest neighbor discrete NLS equation with a 
higher degree of nonlinearity a, which is known Eq] to 
exhibit multistability. 



III. EFFECTS OF ANHARMONICITY: BOUND 

STATE OF DAVYDOV AND BOUSSINESQ 

SOLITONS 

A. System and equations of motion 

Let us consider in this section an excitation energy (or 
excess electron) transport in the anharmonic molecular 
chain. In Davydov's coherent states approximation ]3C| ] 
the Hamiltonian of the system takes on the form 

n 



-1/2 & iM ^ 



n 

2 



Kf) 



• «(a 



'ra+1 



aw(J5; 



'n+l 



/? n ) s 



(25) 



Here, ip n (t) is the excitation wave function of the 
molecule at site n, /3 n (i) is the displacement of the n- 
th molecule from equilibrium, J is the matrix element 
of the excitation transition, the constant \ characterizes 
the exciton - displacement interaction, M is the mass of 
the molecule, and the parameters u> and a characterize 
the elasticity and the anharmonicity of the lattice, re- 
spectively. 

The equations of motion for ip n (t) and (3 n (t) are 



in^Mt) 



M M2Mt) = u(J3 n +l 



-J(lp n+ 1 + l/Vi-1 - 2tp n ) 

/3 n _i - 2/?„) 
x(l-a(/3 n+ i -/3 n _i)) 

X (|^+i| 2 -k/Wi| 2 )- 



(26) 



The Hamiltonian ( pq ) conserves the number of excita- 
tions in the chain. We shall assume that there is only 
one excitation. Thus the normalization condition for the 
wave function ip n (t) is 



IV'nWI 5 



1 . 



(27) 



Using the continuum limit ip n (t) 



-Jkx 



4>(x,t) and 

Pn{t) — ¥ P(^,t), where x — n£ and £ is the lattice con- 
stant, we shall consider the solutions in the form of trav- 
elling waves 



#E,t) = (8A)- 

d/3/dx = -(Jcos(k£)/x) u(<9) , 



where A = 3x 2 /(wJcos(fc£)) and § = {x 
result we obtain the system of equations 



dti 2 



(p — Aip + 2wp = , 
(u - ABu + gu 2 - ip 2 ) = , 



(28) 

vt)/£. As a 

(29) 
(30) 



where the parameters A , B and g are given by the ex- 
pressions 



A = (2 J + HA)/( J cos(fcf)) - 2 
B = 3(v 2 /v 2 - 1) , 
g = (12aJcos(AO)/x . 



(31) 



Here v — £^/ujJM is the sound velocity, v — 
(2J£/Ti) sm(k£) is the soliton velocity, and the dots de- 
note the differentiation with respect to -d . We shall as- 
sume that the effective mass of excitation is positive, i.e. 
J > 0, and consider the carrier wave vector k in the 
interval < k£ < tt/2. 

At the boundary conditions 

w(±oo) = m(±oo) = ii(±oo) = (^(±00) = tp(±oo) = , 



the equations of motion (29) and (|3C|) are similar to that 
of the Henon-Heiles system J41J except for the sign in 
front to ip 2 in Eq. (BO). This similarity has been used 
in our investigations |3q , pq , |39|| of the system — we shall 
review the results in what follows. An alternative ap- 
proach to the problem has been developed recently by 
Zolotaryuk, Spatschek, and Savin |37JJ38|]. 



B. The soliton solutions in the completely integrable 

case 

It is well known that in a general case the Henon-Heiles 
system [|41[ is not completely integrable. However in the 
following three cases |E§] : 

(i)A = 4B,g = l, 
(ii) AA = B , g = 16 , 

(iii) 3 = 6 and arbitrary A and B , 

there exists a second integral of motion, and thus it is a 
Liouville completely integrable system. The conditions 
(i) and (ii) can be satisfied for only one value of the soli- 
ton velocity since the normalization condition ( p7J ) im- 
poses a link between parameters A and B. But we intend 
to investigate how the soliton energy and shape depend 
on the soliton velocity. Therefore we shall consider here 
only the third case. It was shown in Refs. [p5],[36| that in 
this case there exist only three types of soliton solutions 
of Eqs. <M) and (Wh, namely: 



(a) The Boussinesq soliton: 



if = , u——B sech vBd . 
9 



(32) 



It exists at supersonic velocities for any g, and repre- 
sents a lattice compression which moves along the chain 
without changing its shape. 

(b) The Davydov soliton [||: 

ip = 2^/a(A~B) sechVAi} , u = Asech 2 VAtf . (33) 

It exists at subsonic {B < 0) as well as at supersonic 
(B > 0) velocities, but only for g = 6. 

(c) The two-bell shaped soliton, which exists only at 
supersonic (B > 0) velocities and 5 = 6: 



ip = 2VA(A - B)S~ 1 {tf, R) cosh"/B(tf - R) , 

" = ^5(0,*), 

where R is an integration constant and 



(34) 



S(tf, R) = VA cosh(VB(§ - R)) cosh(VAtf) 

- y/Bsmh(VB(d - R)) amh(VM) . (35) 

This solution was found analytically and studied numeri- 
cally in Refs. |p5|p6[ , where it was interpreted as a bound 
state of Davydov and Boussinesq solitons. It was also 
found numerically in Refs. |37],[58| where the authors be- 
lieve, in contrast to our conclusion, that it is a bound 
state of a Davydov soliton with two Boussinesq solitons. 
It is interesting to note that the deformation function 
( p5| ) has the form of the two-soliton solution of the KdV- 
equation Jl9|. So we can conclude that the excitation 
y(#) is similar to that of a quantum particle moving in 
two- well potential — 2u($) (see Figure g). One of the 
wells is created by the lattice soliton, and the second is 
caused by the interaction of the excitation with the lat- 
tice. Part of the time the particle lives in the well that 
was dug by itself, then it tunnels to the well that was cre- 
ated by the lattice soliton, and so on. As a result of such 
a complicated behavior we obtain a wave function y(i?) 
in the form of Eq. ( B4J) (see Figure ||) . It is interesting to 
remark that at the velocity within the range 



i; <«<t'o[l + (A 2 /54) 1/3 ] 1/2 



(36) 



the wave function of the soliton has the one-bell shape 
at any distance R between the wells in the poten- 
tial — 2u('&). On the contrary, if the velocity exceeds 
v a [l + (A/2) 2 / 3 ] 1 / 2 we obtain the two-bell shape at any 
R. However, the ratio of the soliton maximal values in 
this case is proportional to sechVAR at \R\ ^> 1/y/B. 
Thus we may conclude that the excitation (or electron) 
lives mostly in the well that was dug out by itself, tun- 
nelling to the other well with an exponentially small rate. 
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FIG. 6. The form of the two-bell shaped solitons at A = 1 
and 11 = 1.2«o for different values of distance R: R — (solid), 
R — 1 (dashed), and R — 6 (dotted). 

We have calculated (see Figure [7]) the velocity depen- 
dence of the soliton energy for all three types of solitons 
at g = 6. For the energy of Boussinesq soliton ( p2| ) we 
get 



r„.-£(i + «*)*< = 



^Bq 



(37) 



The energy of Davydov soliton (|33|) is given by the ex- 
pression 



Edclv — E e 



where 



Jcos(ki) 



2 
5A 



2 ( l + A v \ A m„L A 
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ex ~ %jp 



cos(k£) 



(38) 



(39) 



is the exciton energy and A is determined by the equation 

VA{A -B)=A. (40) 

Finally the energy of the two-bell shaped soliton turned 
out to be an exact sum of the energies of the Davydov 
and Boussinesq solitons: 



Ell = Eoav + Ebq 



(41) 



It is interesting that this energy does not depend on the 
distance R. It means that at g = 6 the Davydov and 
Boussinesq solitons do not interact. However, as it will 
be shown below they interact in the non-integrable case 



g 7^ 6. Namely, they repel each other at g < 6 and form 
a bound state at g > 6. 



d 2 R 

dr 2 



sinh(VBiJ) 



B cosh^V BR) 



= 0, 



(44) 
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FIG. 7. The energy of the Davydov (Edcv), Boussinesq 
{Esq) and two-bell shaped (En) solitons vs the velocity, for 
A = 10. 



C. Variational approach in the non-integrable case 

In what follows we shall develop a variational approach 
to the investigation of the behavior of the two-bell shaped 
solitons in the anharmonic chain far from the completely 
integrable case. Considering the Hamiltonian (|25|), one 
can see that the Lagrangian of the chain has the form 



n 



-x(/3„+i 

n 



Pn-lH'n 
2 



2a 



{f3 n+l -p n y~— (p n+l -/3 n y 



(42) 



Substituting into it the solution ( |34| ) of Eqs. ( p9[ ) and 
( pfj| ) as a trial function (A and B are supposed to be con- 
stant and R is time dependent), we obtain an effective 
Lagrangian L(R,dtR). If the soliton velocity is close to 
the sound velocity (B <C A) the energy of the two-bell 
shaped soliton is 



E 



ii 



E 



Dav 



E 



Bq 



J cos(M) 



B 3 I' 2 ( dR\ 



dr) 



+Jcos(k£) 



ff-6 

18 



: A 2 / 3 + S{4-sech 2 (v / Bi?)} 



(43) 



The equation of motion in this case is 



where v 2 = (g - 6)VSA/18 and r = y/u/Mt is the 
dimensionless time. 

One can see that at g < 6 there is an intersoliton re- 
pulsion. The distance R — ► oo when r — > oo, and the 
two-bell shaped soliton is dissociated into Davydov and 
Boussinesq solitons. 

On the contrary, at g > 6 there is a potential well in 
Eq. (|43]), and the distance R between solitons is a peri- 
odic function 



2tt 
sinh{VBR) = smh(\/BR ) sin — (r - r ) , (45) 



with the period of oscillations 



2tt 



T = — cosh(VSi? ) 



(46) 



Here Ro is the maximum value of R and To is the initial 
time. 
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FIG. 8. Numerical calculations of the dynamics of two-bell 
shaped solitons for g = 8.5 at different moments of dimension- 
less time: (a) r = 0, (b) r = 1550, (c) r = 3510, (d) r = 5400, 
(c) r = 6460. The boundary conditions are periodic. 

Thus, we can conclude that in highly anharmonic 
molecular chains, there should exist a bound state of 
Davydov and Boussinesq solitons. The dynamical sta- 
bility of this bound state has been tested numerically. 
In Figure @ we present the results of numerical simu- 
lations for g = 8.5. As an initial state we chose an 
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asymmetric two-bell shaped soliton. One can see that as 
time increases, the bells approach each other and their 
heights (as well as the depths of wells) equalize. There- 
after, the distance between the bells increases, and a mir- 
ror reflected asymmetric two-bell shaped soliton appears. 
Thus, the Davydov and Boussinesq solitons form a rather 
stable bound state and demonstrate an oscillating motion 
in accordance with Eq. ( |45| ) . 

In conclusion, we note that for the Davydov vibrational 
soliton motion in the a-helix molecule the parameter 
g ~ 0.8, so that the bound state of supersonic Davydov 
and Boussinesq soliton cannot exist. But for the elec- 
tron motion in the polypeptide molecule we get g > 12. 
Thus we arrive at the conclusion that the bound state 
of the supersonic Davydov electrosoliton and the Boussi- 
nesq soliton can exist in the a-helix biopolymers. 



when the value of the effective (dimensionless) anhar- 
monicity parameter is large enough, a bound state of 
Davydov and Boussinesq solitons can exist. For the 
Davydov vibrational soliton motion in the a-helix pro- 
teins the anharmonicity parameter is too small, so that 
the bound state cannot exist. But for the excess elec- 
tron motion in the polypeptide molecule it is rather large, 
and a bound state of the supersonic Davydov electrosoli- 
ton and the Boussinesq soliton can be formed. Further- 
more, a slightly generalized version of the DNA model 
E3 which takes into account next neighbors transversal 
interactions between Toda chains actually yields strong 
anharmonicity. Hence there is a wide range of values for 
the anharmonicity parameter g occurring in real systems. 
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IV. CONCLUSION 

We have considered two models for energy and charge 
transport and storage in biomolecules. In these mod- 
els we took into account the long-range dispersive (first 
model) and anharmonic (second model) interactions, and 
showed that these interactions are responsible for the ex- 
istence of new types of excitations. 

We have proposed a new nonlocal discrete nonlinear 
Schrodinger (NLS) model for the modelling of the non- 
linear dynamics of the DNA molecule with long-range 
(r~ s and e~@ r ) dispersive interaction between its charged 
groups. We have shown that when the long-range interac- 
tions decay slowly, there is an energy interval where two 
stable stationary states exist at each value of the Hamil- 
tonian H. One of these states is a continuum-like soliton 
and the other one is an intrinsically localized mode. By 
this means a bistability phenomenon came into existence 
in the system. This phenomenon is a result of the com- 
petition of two length scales: the usual scale of the NLS 
model, which is related to the competition between non- 
linearity and dispersion (expressed in terms of the ratio 
N/J), and the radius of the long-range interactions. 

We have shown that a controlled switching between 
narrow (pinned) states and broad (mobile) states is pos- 
sible. Applying a perturbation in the form of a paramet- 
ric kick, we demonstrated that switching occurs beyond 
some well-defined threshold value of the kick strength. 
The particular choice of perturbation is not important 
for the qualitative features of the switching, as long as 
there is a substantial overlap between the perturbation 
and the internal breathing mode. Thus, we believe that 
the mechanism for switching described here can be ap- 
plied for any multistable system where the instability is 
connected with a breathing mode. The switching phe- 
nomenon could be important for controlling energy stor- 
age and transport in DNA molecules. 

The second model is proposed for modelling of anhar- 
monic biological macromolecules. We have shown that 



Yu.G. and S.M. acknowledge support from the 
Ukrainian Fundamental Research Foundation (Grant 
No. 2.4/355). Yu.G. acknowledges also partial financial 
support from SRC QM "Vidhuk". M.J. acknowledges 
financial support from the Swedish Foundation STINT. 



[1] Reiss, C, in M. Peyrard (ed.) Nonlinear Excitations in 
Biomolecules, Springer- Verlag, Berlin, Heidelberg, Les 
Editions de Physique Les Ulis, 1995, p. 29. 
[2] Lozansky et al, in R.H. Sarma (ed.) Stereodynamics of 

Molecular Systems, Pergamon Press, 1979, p. 265. 
[3] Chou, K.C. and Mao, B., Biopolymers, 27 (1988), 1795. 
[4] Prohovsky, E.W., in R.H. Sarma and M.H. Sarma (eds.) 
Biomolecular Stereodynamics IV, Adenine Guilderland, 
New York, 1986, p. 21. 
[5] Georghiou et al, Biophysical J., 70 (1996), 1909. 
[6] Holbrook, SR. and Kim, S.H., J. Mol. Biol., 173 (1984), 

361. 
[7] Gaididei, Yu.B. and Serikov, A.A., Theor. and Math. 

Phys., 27 (1976), 457. 
[8] Papanicolaou, G.C., J. Appl. Math., 21 (1971), 13. 
[9] Wadati, M., J. Phys. Soc. Jpn., 38 (1976), 673. 
[10] Sievers, A.J. and Takeno, S., Phys. Rev. Lett., 61 (1988), 

970. 
[11] MacKay, R.S. and Aubry, S., Nonlinearity, 7 (1994), 

1623. 
[12] Pouget, J. et al, Phys. Rev. B, 47 (1993), 14866. 
[13] Dauxois, T. and Peyrard, M., Phys. Rev. Lett., 70 (1993), 

3935. 
[14] Muto, V. et al, Phys. Rev. A, 42 (1990), 7452. 
[15] Christiansen, PL. et al, Phys. Rev. B, 55 (1997), 5759. 
[16] Gaeta, G et al, Riv. N. China., 17 (4) (1994), 1. 
[17] Peyrard, M. and Bishop, A., Phys. Rev. Lett., 62 (1989), 

2755. 
[18] Techera, M. et al, Phys. Rev. A, 40 (1989), 6636. 
[19] Dauxois, T. et al, Phys. Rev. E, 47 (1993), 684. 



11 



[20] Dauxois, T. et al, Phys. Rev. E, 47 (1993), R44. 

[21] Dahlborg, U. and Rupprecht, A., Biopolymers, 10 (1971), 

849. 
[22] Corongiu, G. and Clementi, E., Biopolymers, 20 (1981), 

551. 
[23] Braun, O.M. et al, Phys. Rev. B, 41 (1990), 7118. 
[24] Woafo, P. et al, J. Phys. Condens. Matter, 5 (1993), L123. 
[25] Vazquez, L. et al, Phys. Lett. A, 189 (1994), 454. 
[26] Alfimov, G.L et al, Chaos, 3 (1993), 405. 
[27] Gaididei, Yu.B. et al, Phys. Lett. A 222 (1996), 152; 

Gaididci, Yu.B. et al, Phys. Scr. T67 (1996), 151. 
[28] Gaididei, Yu.B. et al, Phys. Rev. Lett., 75 (1995), 2240. 
[29] Gaididei, Yu.B. et al, Phys. Rev. E, 55 (1997), 6141. 
[30] Davydov, A.S.: Solitons in molecular systems, D. Reidel 

Pub. Co., Dordrecht, Boston, Lancaster, 1987. 
[31] Christiansen, P.L. and Scott, A.C. (eds.): Davydov's 

soliton revisited, NATO ASI Series B: Physics Vol. 243, 

Plenum Press, 1990. 
[32] Christiansen, P.L. et al, Chaos, Solitons and Fractals, 2 

(1992), 45. 
[33] Davydov, A.S. and Zolotariuk, A.V. Phys. Stat. Sol. (b), 

115 (1983), 115. 
[34] Davydov, A.S. and Zolotariuk, A.V., Phys. Scr., 30 

(1984), 426. 



[35] 
[36] 
[37] 
[38] 

[39] 



[40] 

[41] 
[42] 

[43] 

[44] 
[45] 
[46] 

[47] 

[48] 



[49] 



Christiansen, P.L. et al, Phys. Lett. A, 166 (1992), 129. 

Gaididei, Yu.B. et al, Phys. Scr., 51 (1995), 289. 

Zolotaryuk, A.V. et al, Europhys. Lett., 31 (1995), 531. 

Zolotaryuk, A.V. et al, Phys. Rev. B, 54 (1996), 266. 

Mingaleev, S.F., Effects of Nonlocality and Anharmonic- 

ity in Nonlinear transport of energy and charge, Ph.D. 

thesis, (Kiev, 1997). 

Perez, P. and Theodorakopoulos, N., Phys. Lett. A, 124 

(1987), 267. 

Henon, M. and Heiles, C, Astrophysics, 63 (1964), 73. 

Baker, G.A. Jr, Phys. Rev., 122 (1961), 1477. 

Kac, A.M. and Helfand, B.C., J. Math. Phys., 4 (1972), 

1078. 

Johansson, M. et al, Phys. Rev. E, 57 (1998), 4739. 

Laedke, E.W. et al, Phys. Rev. Lett., 73 (1994), 1055. 

Malomed, B. and Weinstein, M.I., Phys. Lett. A, 220 

(1996), 91. 

Johansson, M. et al, Physica D, 119 (1998), 115. 

Lichtenberg, A.J. and Lieberman, M.A.: Regular and 

Stochastic motion, Springer- Verlag, Berlin, Heidelberg, 

New York, 1983. 

Lamb, G.L. Jr.: Elements of Soliton Theory, John Wiley 

and Sons, New York, 1980. 



12 



